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ABSTRACT 


Magnetic reconnection is an important physical process in various explosive 
phenomena in the universe. In the previous studies, it was found that fast re¬ 
connection takes place when the thickness of a current sheet becomes on the 
order of a microscopic length such as the ion larmor radius or the ion inertial 
length. In this study, we investigated the pinching process of a current sheet by 
the Lorentz force in a \ow-f3 plasma using one-dimensional magnetohydrodynam¬ 
ics (MHD) simulations. It is known that there is an exact self-similar solution 
for this problem that neglects gas pressure. We compared the non-linear MHD 
dynamics with the analytic self-similar solution. From the MHD simulations, we 
found that with the gas pressure included the implosion process deviates from 
the analytic self-similar solution as f —)■ to, where to is the explosion time when 
the thickness of a current sheet of the analytic solution becomes 0. We also found 
a pair of MHD fast-mode shocks are generated and propagate after the formation 
of the pinched current sheet as t —?■ to- On the basis of the Rankine-Hugoniot 
relations, we derived the scaling law of the physical quantities with respect to 
the initial plasma beta in the pinched current sheet. Our study could help us to 
estimate the physical quantities in the pinched current sheet formed in a low-/3 
plasma. 


Subject headings: Magnetohydrodynamics (MHD) - Stars:sun - Stars: magnetic field 
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Introduction 


It has been known that magnetic reconnection plays important roles in ma ny energetic 


explosions like solar flares (IZweibel fc Yamada 


2009 


Shibata fc Magara 


20111). The 


Masada et al. 

2010; 

Meng et al. 

— ^ 

r—1 

o 

CM 


Maehara et al. 

2012 

) and soft gamma repeaters ( 

Lvutikov 

2006 


2014J) . In a classical standard model of solar flares, the 


magnetic energy is rapi dly converted to the thermal anc 


magnetic reconnection fICarmichael 


1976 


Shibata et al 


1995 


1964: 


Sturrock 


Shibata fc Magara 


201 


much shorter current sheet width than its length (lYamada et ah 


1966 


dnetic energies 


Hiravama 


1974 


of plasmas through 


Kopp fc Pneuman 


) ■ Such rapid recon nection requires a 


2010l ). Therefore, how thin 


current sheets are established is a central issue for understanding the origin of explosive 
phenomena. 


Magnetohydrodynamics (MHD) proc esses to form a t 
theoretically studied by many authors. In 


lin current sheet have been 


Forbes & Priest 

(1995 

) and 

Lin & Forbes 

(20001. 


they numerically and analytically investigated the formation of a current sheet below an 
erupting flux rope during a solar flare. There are many observational supports of the 
formation and the evol ution of a curre nt sheet between a cusp-shaped flaring loop and 


an erupting flux rope (ILiu et al. 


studied flForbes fc 


McLaughlin et al. 


201011 . The formation of a current sheet between an 
he solar corona has been nu merically 


a pre-existing magnetic fie 

d in t 

(’riest 

19§4; 

Shibata et al. 

1992; 


Yokoyama fc Shibata 


19961 1. In 


(1200911 . the collapse of a null point to a current sheet triggered by MHD 


fast-mode waves was investigated. 

A long and thin current sheet can form many fine scale current sheets inside it through 
some MHD instabilities. Such a current sheet can be frag mented to form magnetic islands 


or plasmoids through the tearing or plasmoid instabilities fjFnrth et al 


1963 


Lonreiro et al. 
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20071) . The plasmoids are pulled against each other by a Lorentz force. When plasmoids 
collide with each other, new current sheets will be formed between t hem. This process c an 


operate at multi-spatial scale, depenc 


Shibata fc Tanuma 


2001 


Bartaetah 


ingon the Lundquist number (ITajima et ah 


1987 : 


20111) . Observational suppor t from the dire c t imag ing 


observation of the plasmoid formation in a solar flare was given by 


Takasao et al. 


(120121) for 


the hrst time. They also found that the plasma blobs (possibly plasmoids) in the current 
sheet collided or merged with each other before they were ejected from the current sheet. 
The pinching process of a current sheet by a Lorentz force like the last two processes above 
is the main subject of this paper. 

In this study, we investigated the implosion process in which current sheets are pinched 
by the magnetic pressure gradient. In a previous study, the self-similar solution in the 


implosion proce ss was analytically derived under the assumption that 


can 


1987h . In 


3e ig nored fllmshennik fc Svrovatskii 


19671: 


Bulanov fc Olshanetskii 


he gas pressure 


19841: 


Taiima et al 


ForbesI (119821) . they investigated the implosion process including the gas pressure. 


and found the generation of two MHD fast-mode shocks as a result of a signihcant increase 
in the gas pressure at the center of the current sheet. However, in their simulation, the 
current sheet was incompletely pinched because of the low spatial resolution. To study the 
pinching process of the thin current sheet formed behind the shocks in detail, we performed 
a series of numerical simulations with a sufficient spatial resolution. 


In §2, the self-similar solution proposed by 


Taiima et al. 


fjl987() will be introduced. In 


§3, we will show the basic equations and initial conditions of the numerical simulations. 
In §4, we will show the results of the numerical simulations and the power law scalings of 
the physical quantities in the pinched current sheets. In §5, we will analytically derive the 
power law relation shown in §4. Finally, we will summarize our conclusion in §5. 
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Self-Similar Solution 


The situation that we consider is similar to that of iTaiima et al.l fjl987|) where a 
current sheet is pinched by a Lorentz force (see Figure [I]). They derived a self-similar 
solution using the two-fluid equations and Maxwell’s equations. Here we derive the same 
self-similar solution from ID MHD equations to review the basic physics and assumptions 
of the solution. We start from the basic ID MHD equations: 


d 9 ^ 

(pv.) . 

dt ^ dx^ 

■ - ——R 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Now we introduce dimensionless physical quantities as below: 


X* = x/L, 
t* = t/xA, 
p*{t*) = p(t)/po, 
p*{x\t*) =p{x,t)/po, 
B*{x*,t*) = By{x,t)/Bo, 
v*{x*,t*) = v^{x,t)/VA, 


where L is the initial current sheet thickness, Va is the initial Alfven speed, ta = L/Va is 
an Alfven timescale and the initial values of the physical quantities are denoted by 0. The 
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dimensionless equations can be written as 


_d 

dt 


d 


:P = 


dx* 


{p*v*), 


d d d 


{p*v*^) 


A. 

dx* 


Bf 


—B* = 
dt* ^ 


d 


dx* 




(5) 

( 6 ) 

(7) 


where Po = Po/(-Bo/2) is the initial plasma beta. 

We introduce a scale factor a(t) as follows: 

* h(r) , 

V = ——rX , 

a{t*) 


( 8 ) 


where a dot represents the time derivative. An Ansatz is that the velocity is linear in x*. 
For simplicity, we assume that p*, the density in the current sheet, is spatial uniform. From 
equation ([5]) and (0), we obtain 

1 


P* = 


a(t*)' 


(9) 


We also assume that B* is linear in x*] B*{x*,t*) = B*{t*)x*. Then from equation 
and ([HI), we obtain 

B* = ^ 

“ a{t*f 

We take the following form of the scaling factor. 


a{t*) oc 


to -1 
to 


= 1 


* / ’ 


( 10 ) 


( 11 ) 


where to is the explosion time when the thickness of the current sheet goes to 0. Neglecting 
the gas pressure term in ([6]), we obtain the index k and explosion time to from equation (l6|) 
and flTT]) : 


a(t*) oc 1 - - 


f* — 
an — 


2/3 


( 12 ) 


^/2 


( 13 ) 
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From equation ffT^ . we obtain the following expressions: 


v*{x*, t*) oc 



(14) 



2/3 


(15) 



4/3 


oc 1 


(16) 


Note t hat these scaling laws are the same as those original obtained bv llmshennik fc Svrovatskii 


dlMZl). 


If the plasma adiabatically evolves, we obtain the gas pressure in the following form, 

- 10/9 



( 17 ) 


This indicates that the gas pressure at the center should explode as f ^ to- Therefore the 
pinching by a Lorentz force must stop, which cannot be described by the linear theory. 


To study the evolution near and after the explosion time, we numerically investigated the 
implosion of a current sheet using ID MHD simulations. 

3. Basic Equations and Initial Conditions of Numerical Simulations 

All the physical quantities are functions of x and t, and a magnetic held is considered 
only in the y-direction {By). The basic equations are as follows: 



(18) 


(19) 


( 20 ) 


( 21 ) 








where the total energy density e is dehned as 


P , pVx^ , By^ 

e =--- 

7-1 2 Stt 


( 22 ) 


We adopted the self-similar solution of 


Tajima et al. 


fjl987l) as the initial condition of 


the numerical simulations. The initial condition is not an MHD equilibrium, and after the 
simulation’s start, the current sheet is pinched by a Lorentz force as shown in Figure [TJ 
The calculation domain is in the range of 0 < x/L < 2, where L is the initial width of the 
current sheet. An anti-parallel magnetic held is given by 


By{x) 


= Untanh ( — 


X 


(23) 


where Bq is dehned as | By{:kL) \. For simplicity, we set the initial gas pressure p and the 
initial electron density p to be spatially uniform. The initial gas pressure is normalized by 
the initial magnetic pressure Bq^/Sti. The initial velocity held is given by 


Vxix) = uotanh 


X 


(24) 


In our numerical s imulati o ns, t 


larger than 512 in 


le total grid number is hxed to 40000, which is adequately 


ForbesI fjl982l ). The numerical scheme we adoptee 


scheme, which is a fully shock-capturing scheme flMivoshi fc Kusano 


is based on the HLLD 


200, 5h 


Results of Numerical Simulations 


4.1. Time Evolution of a Current Sheet in Nnmerical Simulation 


As a typical case, here we show the time evolution of the c ase with dn = 


Figure [2]). During t < to, the self-similar evolution predicted in 


Tajima et al. 


-2 


see 


fll987h is 


conhrmed with plasma gas pressure (see Figure [3]). At t to, the implosion leads to the 
formation of a shock (MHD fast-mode shock) as a result of a signiheant increase in the gas 
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pressure at the center of the current sheet. As shown in Figure [2], the compressed plasma 
is continuously left behind the shock as the shock propagates outward. Note that p/p^ is 
almost kept constant at the center of the current sheet, which indicates that the physical 
quantities there adiabatically evolve. To investigate the implosion process with the gas 
pressure, we performed the numerical simulations with various values of initial gas pressure 
and the spatial constant value of initial velocity. As a result of numerical simulations, we 
found the power law behavior of physical quantities in the pinched current sheet behind the 
shock. As shown in next section, we were able to analytically derive the scaling law with 
respect to Pq. 

4.2. /^o Dependence of Physical Quantities in Current Sheet 

We investigate the (3q dependence of the physical quantities in the current sheet. The 
range of /So is —3 < log/So < — 1- The /So dependence of the density, pressure, and magnetic 
held strength just behind the shock is shown in Figure 01 Figure 0] implies the following 
power law relations: 


pIpo «. I3„ 
p/p„ oc 


(26) 

(26) 

(27) 


By/Byo ~ const. 


We also investigated the /So dependence of the thickness of a current sheet, Amin- We 
dehne Lrain as follows: 



(28) 
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where and Jmax(^) are respectively 

^max = max{By{x,t )), 

Jmax = iaax{J^{x,t)). 

. As a result of the numerical simulations, we found a scaling law of Lj^in (see Figure [5]) as 
follows: 


Lmin CX (29) 

To study the formation of the pinched current sheet in detail, we tracked selected 
Lagrangian particles that are initially located in the initial current sheet. Figure [6] shows 
the trajectories of the Lagrangian particles on the time-distance diagrams of na;,logioP and 
By. From the Figure El we found that the pinched current sheet is formed from a part of 
the initial current sheet, not from the whole of it. 


5. The Analytical Discussion about the /3o Dependence of Physical Quantities 


In this section, we aim to derive the dependence of the physical quantities in the 


pinched current sheet. As shown in §4, we conhrmed the pinching p rocess 


force and the propagation of the MHD fast-mode shocks, as shown in 


jy the Lorentz 


Forbes 


f 1982 1. Behind 


the fast-mode shock, the physical quantities are determined by the Rankine-Hugoniot 
equations: 


Vx2 

Vxl 

By2 

Ai 

El 

Pi 


1 

X’ 







+ 


1 


+ 1) 


(30) 

(31) 


(32) 
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where the physical quantities in front of and behind the shock are respectively denoted by 
1 and 2, Mi = Vxi/csi is the Mach number, and X = P 2 /Pi is the density ratio, which is the 
positive solution of 


2 (2 - 7) + 2/3i + (7 - 1) /SiMi + 2 


7X 


-7 (7 + 1) l 3 iMi = 0. 


(33) 


From equation (13T1) and (l33|) . the compression ratio is limited in the range of 


B 2 7 + 1 
1 < ^ < 


(34) 


Bi 7 - 1 ’ 

where the upper limit is 4 for 7 = 5/3. Therefore the maximum value of the strength 
of the magnetic held should weakly depend on /^o (i?max/5o ~ const with respect to /do, 
equation fl27j) ). The current sheet, which is pinched by the shocks, is in the MHD equilibrium 
(i.e. Ua; = 0). This means that the gas pressure at the center of the current sheet, Pmax, 
should be comparable to the magnetic pressure outside of the current sheet: 


By, 


B. 


Pn 


''max -^0 

r\j - r\j 


Stt 


Sir ’ 


(35) 


which gives the relation fl26|) 


Pmax 

Po 


P<r\ 


(36) 


where po is the initial gas pressure. From equation fl36ll and the adiabatic condition, the 
density at the center (which takes the maximum value, equation fl25lF can be written by 


Pmax / Pmax 

- r\j 


Po V Po 


1/7 


(37) 


where po is the initial plasma density and 7 = 3/5 is the adiabatic index. 

To derive the /do dependence of Tmin, the relation fl29l) . we dehned i?i, the magnetic 
held strength of the maximum value in the pinched current sheet. As shown in Lagrangian 
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particle trajectories, note that the pinched current sheet is formed from not the whole of the 
initial current but a part of it. We also defined the length of the part in the initial current 
sheet as Lq and the maximum value of the magnetic held strength in the part as Bq. In our 
numerical simulations, we conhrmed Lq Lq by the Lagrangian particle motions. Near the 
center of the initial current sheet, the magnetic held strength linearly depend on x. 


Bo' 


We considered the conservation law of a magnetic hux. 


(38) 


B*Ll = BiL^ 


= BoL 


min* 


From equation fl38l) and ([39]), we obtained 



Finally, from the mass conservation law. 


Lq 

Lmin 


Pmax 

Po 



From equation fHOl) and (HTl) . we derived 


o6/5 
~ PO 5 

which is consistent with equation fl2^ . 


(39) 


(40) 


(41) 


(42) 


6. Discussion 


In this paper, we studied the implosion process of current sheets in a low-/3 plasma 
usi ng one-dimen s ional ideal MHD simulations. We conhrmed that the self-similar solution 


by 


Taiima et al. 


fll987[) holds before the explosion time. The plasma adiabatically evolves 
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till the explosion time. At t ~ to; the behavior of the current sheet deviates from the 
analytical self-similar solution and the MHD fast-mode shocks are formed because of the 
increase in the gas pressure at the center of the current sheet. After a pair of MHD 
fast-mode shocks propagate away from the current sheet, the pinched current sheet is 
formed between a pair of the shocks. We studied the /3o dependence of the maximum values 
of the physical quantities inside the current sheet and the minimum thickness of the current 
sheet (equation fl36l) . flTTI) and fH2l) ). 


Contrast to our study, the rj dep endence of the physical qua ntities and of the 


reconnection rate were investigated by 


McClvmont fc CraigI f 19961 ) in MHD simulations 


similar to ours, where 7] is the magnetic diffusivity. They found that the effect of the gas 
pressure does not signihcantly change the dependence of the thickness of the current sheet 
on T], which means that our scaling laws can be easily extended to the cases with the 
resistivity. 


In 


McLaughlin et ah 


(120091) . the collapse of a null point to a current sheet triggered 


by MHD fast-mode waves was investigated. In their simulations, successive current sheet 
pinching was observed. The pinching is done by the gas pressure of the heated plasma in the 
outflow region, while the pinching in our study is done by a Lorentz force. To extend this 
study, we are currently performing two-dimensional MHD simulations of the coalescence 
process of plasmoids and studying the pinching process of the current sheet formed between 
plasmoids in detail. 


We analyzed the /So dependence of the fast-mode Mach number and the compression 
ratio, p 2 /pi{see Figure [T]). Figure [7] shows that the fast-mode Mach number. Mi, weakly 
depend on (3q. We also found that the compression ratio of the density and magnetic held. 
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X, weakly depend on (3q. Therefore, from eqnation ([32]), we obtain 


■oSio I ^ I oc A ~ /3o- 


(43) 


This /So dependence of the compression ratio, P 2 /P 1 , is similar to eqnation fl45D with /^q < 10 ^ 
11 recent stndies, it is indicated that particles can be accelerat ed in contracting plasmoid 


Drake et ah 


20061) and in plasmoids crossing fast-mode shocks IjNishiznka fc Shibata 


2OI3I1 . 


Fignre [Hl(a) shows an overview of the implosion process of the current sheet between 
plasmoids and the formation of the MHD fast-mode shocks. In our study, we found that 
fast-mode shocks can be formed during the implosion process in low-/9 plasmas. Applying 
the shock formation in the implosion process to a coalescence process, we conjecture that 
particles in the plasmoids are efficiently accelerated at shocks generated by the coalescence 
of plasmoids through a Fermi acceleration process (see Figure [8] (b)). 

We investigated the parameter region where our scaling laws are valid. Since we assume 
the framework of non-relativistic ideal MHD, (1) the drift velocity of the electrons needs to 
be sufficiently smaller than the light speed and (2) the thickness of the current sheet needs 
to be larger than the ion skin depth. The electron drift velocity can be written as follows: 


c ^ „ c B 

= q—V X B 

47r 47r L 

where L is the thickness of current sheet. From equation 
electrons, is given as 

Wrift,e f B 

- r\j -^ 

c 47re riL 

Applying the scaling laws of the relation (l36l) . (ITT]) and (IT2|) to equation 

Wrift.e 


(44) 

, fdrift.e, the drift velocity of 


(45) 


we obtain 


10 


-4 


/So fBo\ fnoY^ [Lo 




Bi 


ng 


(46) 

















15 


where /1_2 = 10 Bi = lOG, ng = lO^cm ^ and Lg = 10®cm. Since we study non-relativistic 
MHD processes, ndrift,e/c < 1 is required. This leads to 


logio 



> - logio 
+ logio 




(47) 


Since the thickness of the current sheet must be larger than the ion skin depth, we 
similarly obtain 


d 


1 

c 



10 “^ 




(48) 


where Aj is the ion skin depth. We can transform equation fl48l) as 


logio 



>2 logio 




-4. 


(49) 


In Figure [9l we show the parameter region using inequality (1471) and with the value 
in the solar corona, /Ig = 10“^ and Bq = lOG. It is found that our scaling laws are applicable 
in the solar corona. These scaling laws will help us to estimate physical quantities in a 
current sheet formed in a low-/9 plasma. Describing the evolution of the pinching process in 
an extremely low-/? plasma like the atmosphere of the neutron stars will be our future work. 
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acknowledges support by the Research Fellowship of the Japan Society for the Promotion 
of Science (JSPS). This work was supported by the Grant-in-Aids from the Ministry of 
Education, Gulture, Sports, Science and Technology of Japan (No. 25287039). 







- 16 - 

REFERENCES 

Barta, M., Buchner, J., Karlicky, M., & Skala, J. 2011, ApJ, 737, 24 

Bulanov, S. V., & Olshanetskii, M. A. 1984, Physics Letters A, 100, 35 

Carmichael, H. 1964, NASA Special Publication, 50, 451 

Drake, J. F., Swisdak, M., Che, H., & Shay, M. A. 2006, Nature, 443, 553 

Forbes, T. G. 1982, Journal of Plasma Physics, 27, 491 

Forbes, T. G., & Priest, E. R. 1984, Sol. Phys., 94, 315 

Forbes, T. G., & Priest, E. R. 1995, ApJ, 446, 377 

Furth, H. P., Killeen, J., & Rosenbluth, M. N. 1963, Physics of Fluids, 6, 459 
Hirayama, T. 1974, Sol. Phys., 34, 323 

Imshennik, V. S., & Syrovatskii, S. I. 1967, Soviet Journal of Experimental and Theoretical 
Physics, 25, 656 

Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys., 50, 85 
Lin, J., & Forbes, T. G. 2000, J. Geophys. Res., 105, 2375 
Liu, R., Lee, J., Wang, T., et ah 2010, ApJ, 723, L28 

Loureiro, N. F., Schekochihin, A. A., &: Cowley, S. C. 2007, Physics of Plasmas, 14, 100703 
Lyutikov, M. 2006, MNRAS, 367, 1594 

Maehara, H., Shibayama, T., Notsu, S., et al. 2012, Nature, 485, 478 
Masada, Y., Nagataki, S., Shibata, K., & Terasawa, T. 2010, PASJ, 62, 1093 


17 


McClymont, A. N., & Craig, 1. J. D. 1996, ApJ, 466, 487 

McLaughlin, J. A., De Moortel, L, Hood, A. W., & Brady, C. S. 2009, A&A, 493, 227 

Meng, Y., Lin, J., Zhang, L., et ah 2014, ApJ, 785, 62 

Miyoshi, T., & Kusano, K. 2005, Journal of Computational Physics, 208, 315 

Nishizuka, N., & Shibata, K. 2013, Physical Review Letters, 110, 051101 

Shibata, K., & Magara, T. 2011, Living Reviews in Solar Physics, 8, 6 

Shibata, K., Masuda, S., Shimojo, M., et ah 1995, ApJ, 451, L83 

Shibata, K., Nozawa, S., & Matsumoto, R. 1992, PASJ, 44, 265 

Shibata, K., & Tanuma, S. 2001, Earth, Planets, and Space, 53, 473 

Sturrock, P. A. 1966, Nature, 211, 695 

Tajima, T., Sakai, J., Nakajima, H., et al. 1987, ApJ, 321, 1031 
Takasao, S., Asai, A., Isobe, H., & Shibata, K. 2012, ApJ, 745, L6 
Yamada, M., Kulsrud, R., & Ji, H. 2010, Reviews of Modern Physics, 82, 603 
Yokoyama, T., & Shibata, K. 1996, PASJ, 48, 353 
Zweibel, E. G., & Yamada, M. 2009, ARA&A, 47, 291 


This manuscript was prepared with the AAS IATeX macros v5.2. 



18 


(a) 
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y' 
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4 
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X 

N. 

> 

f 
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^ is dorr 
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/ 
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is dominant 

(b) 



Fig. 1.— (a) A schematic picture of the initial condition of our simulations. The vertical 
solid arrows indicate magnetic held lines and the horizontal arrows indicate the plasma how. 
(b) The initial condition of a typical case of our simulations. The solid line represents the 
initial distribution of By and the dashed line represents the initial distribution of v^- Va 
represents the Alfven speed. 
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log ,0 (p/po) 




Fig. 2.— The time evolution of the implosion process of p and By with /3o = 10 The time 
interval and the last time are respectively At = 0.9 and tend = 1-5 in the unit of the Alfven 


time. 
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1 



2 


3 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

logio(l-t/to) 


Fig. 3.— Time evolution of p, p, By and near the center of the current sheet with (3^ = 10 


-2 


The initial values are denoted by subscript 0. In the left column, the diamonds represent the 
nu merical result s and the solid lines represents the analytical self-similar solution proposed 


m 


Taiima et ah 


(jl987|l . In the right column, he diamonds represent the numerical results 


and the solid lines represents the fitting lines by the least squares method. The indexes of the 


result of fitting are quite similar to them of the self-similar solution proposed in 
( 19871 1. 


Taiima et ah 
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-3.0 -2.5 -2.0 -1.5 -1.0 

logio(P o) 


Fig. 4.— The /3o dependence of Pmax,Pmax and i?max- The diamonds represent the result of 
the numerical simulations and the solid lines represent the fitting lines by the least squares 


method. 
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Fig. 5.— (a)The /So dependence of the maximum current density Jmax(top) and of the 
minimum thickness of the current sheet Lmin(bottom). The diamonds represent the result of 
the numerical simulations and the solid lines represent the htting lines by the least squares 
method, (b) The current density distribution with /So = 0.01 and = 1.8, where tA is the 
Alfven time. The vertical line indicates the current sheet thickness, Lqs = -Bmax/'^max- 
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Fig. 6.— Time-distance diagrams of Vx (top), log^gp(niiddle) and i?j/(bottom). The color 
contour denotes the values of the quantities. The solid lines represent the trajectories of the 
selected Lagrangian particles. The vertical dashed lines represent the current sheet thickness, 
Lcs = -Bmax/</max- The left column shows region of 0 < x < 1 and 0 < f < 2. The right 
column shows the limited region of —4 < log^gX < —0.5 and 1.0 <t< 1.4. 


































Mach number 
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Fig. 7.— The /Sq dependence of the fast-mode Mach number and P 2 /P 1 , where the gas 
pressures behind and in front of the shock are respectively denoted by the subscript 1 and 


2. 
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(a) 


to<ti<t2 


o o 




MHD fast-mode shocks 



Fig. 8.— (a)A schematic picture of the time evolution of the pinching process and the 
formation of MHD fast-mode shocks. (b)A schematic picture of our conjecture about an 
actual high-Lundquist number reconnection process, where particles could be accelerated at 
fast-mode shocks generated by the coalescence of plasmoids. 
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-4 -2 0 2 4 -4 -2 0 2 4 -4 -2 0 2 4 

logioK/ng) logio(no/ng) logjoK/ng) 

Fig. 9.— The parameter regions where both inequality (IT7|) and fHOD are valid (i.e. the non- 
relativistic MHD approximation is valid). The cases with /^o = 10“^(left), /Sq = 10“^(middle) 
and /3o = 10“®(right) are shown , where rig = 10®cm“^ and Lg = 10®cm. The shaded zones 
represent the parameter regions where our scaling laws do not hold. 



















































